A generalized Damour-Navier-Stokes equation applied to trapping horizons 
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An identity is derived from Einstein equation for any hypersurface Ti which can be foliated by 
spacelike two-dimensional surfaces. In the case where the hypersurface is null, this identity coin- 
cides with the two-dimensional Navier-Stokes-like equation obtained by Damour in the membrane 
approach to a black hole event horizon. In the case where TC is spacelike or null and the 2-surfaces 
are marginally trapped, this identity applies to Hayward's trapping horizons and to the related 
dynamical horizons recently introduced by Ashtekar and Krishnan. The identity involves a normal 
fundamental form (normal connection 1-form) of the 2-surface, which can be viewed as a general- 
ization to non-null hypersurfaces of the Hajicek 1-form used by Damour. This 1-form is also used 
to define the angular momentum of the horizon. The generalized Damour-Navier-Stokes equation 
leads then to a simple evolution equation for the angular momentum. 

PACS numbers: 04.70.Bw,04.20.-q,04.30.Db,04.25.Dm 



I. INTRODUCTION 

The concept of black hole shear viscosity has been 
introduced by Hawking and Hartle [1-3] when study- 
ing the response of the event horizon to external per- 
turbations. It was then greatly enhanced by Damour 
[4-6] who showed that a 2-dimensional spacelike section 
of the event horizon can be considered as a fluid bub- 
ble endowed with some mechanical and electromagnetic 
properties. Moreover, he explicitly derived from Einstein 
equation a Navier-Stokes equation, involving both shear 
and bulk viscosities, for an effective momentum density 
of the "fluid" constituting the bubble. This "fluid bub- 
ble" point of view lead to the development of the so-called 
membrane paradigm for black hole astrophysics [7, 8] . An 
example of recent work using the concept of black hole 
viscosity is the study of tidal interaction in binary-black- 
hole inspiral [9]. 

The membrane approach was related to the event hori- 
zon of the black hole (or to the associated "stretched 
horizon" [7, 8]). However the event horizon is an ex- 
tremely global and teleological concept, which requires 
the knowledge of the complete spacetime, including the 
full future of any Cauchy surface, to be located. It can by 
no means be determined from local measurements (see 
§ 2.2.2 of Ref. [10] for an interesting example of some 
event horizon which appears in a flat region of space- 
time) . This makes the event horizon a not very practical 
representation of black holes for studies beyond the sta- 
tionary regime in numerical relativity or quantum grav- 
ity. For this reason, local characterizations of black holes 
have been introduced in the last decade (see [10-12] for 
a review). Although some local concepts, like Hawking's 
apparent horizon [13], have appeared well before, the lo- 
cal approach has really started with Hayward's introduc- 
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tion of trapping horizons (or more precisely future outer 
trapping horizons) [14]. Whereas apparent horizons are 
2-dimensional surfaces (associated with some spacelike 
slicing of spacetime), trapping horizons are 3-dimensional 
submanifolds of spacetime (hypersurfaces), as event hori- 
zons. Basically a trapping horizon is a world tube made 
of marginally trapped 2-surfaces; Hayward has studied 
the dynamics of these objects on their own, without mak- 
ing any reference to any slicing of spacetime by spacelike 
Cauchy surfaces. More recently, Ashtekar and Krishnan 
[10, 15, 16] have introduced the related concept of dy- 
namical horizons and have established the "first law" of 
black hole thermodynamics for them. This first law has 
been extended to trapping horizons [17, 18]. 

A natural question which arises is then: can the "fluid 
bubble" approach to event horizons be extended to these 
local characterizations of black holes ? In particular, 
can one obtain an analog of Damour's Navier-Stokes 
equation ? Although some viscosity aspects are already 
present (in form of dissipations terms) in the first law 
of dynamical horizons established by Ashtekar, Krishnan 
and Hayward, it does not seem a priori obvious to get an 
equivalent of the Navier-Stokes equation. In particular, 
Damour's derivation relied heavily on the null structure 
of the event horizon, whereas the future outer trapping 
horizons are generically spacelike in dynamical situations, 
being null only in stationary states. 

We will show here that it is indeed possible to get a 
Navier-Stokes-like equation, provided one introduces the 
correct geometrical objects. Actually, the Navier-Stokes- 
like equation derived hereafter is quite general: it applies 
not only to trapping or dynamical horizons, but to any 
hypersurface which can be foliated by a smooth family of 
spacelike 2-surfaces. In this respect, the recent demon- 
stration [19] of the uniqueness of the foliation of a given 
dynamical horizon by marginally trapped surfaces is pro- 
viding some motivation for the present work. Moreover, 
some uniqueness theorems about dynamical and trapping 
horizons have been recently established [19, 20], confer- 



ring to these objects a more solid physical status. In 
addition, the recent study [21] has provided some deep 
insights about their behavior in various scenarios of grav- 
itational collapse. 

The plan of the article is as follows. In Sec. II, we set 
the basic framework of our study, namely a hypersurface 
foliated by spacelike 2-surfaces. In Sec. Ill we review 
standard results about the extrinsic geometry of a single 
spacelike 2-surface. At this stage we introduce the geo- 
metrical object that shall play the role of a momentum 
density in the Navier-Stokes equation, namely a normal 
fundamental form of the 2-surface. In Sec. IV are intro- 
duced other geometrical objects, defined by the foliation 
as a whole and not by a single 2-surface. Then we have all 
the tools to derive the generalized Damour-Navier-Stokes 
equation in Sec. V. An application to a general law of 
angular momentum balance is given in Sec. VI. Finally 
Sec. VII contains the concluding remarks. 
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FIG. 1: Foliation of a hypersurface TL by a family (Stjtev, of 
spacelike 2-surfaces, and the associated evolution vector h. 



II. FOLIATION OF A HYPERSURFACE BY 
SPACELIKE 2-SURFACES 

A. General set-up 

We consider a spacetime (A4,g), i.e. a smooth mani- 
fold M. of dimension 4 endowed with a Lorentzian met- 
ric g, of signature (—,+,+,+). We assume that M. is 
time-orientable. Let TL be a hypersurface of M. which is 
foliated by a family (St)teR of 2-dimensional surfaces St 
labeled by the real parameter t. By foliation, it is meant 
that TL = UteK &t an d that for each point p e TL, there 
is only one St going through p (see Fig. 1). Then, given 
a coordinate system x a = (x 2 ,x 3 ) on each St, (t,x 2 ,x 3 ) 
constitutes a coordinate system on TL. 1 We assume 
that all surfaces St are spacelike and closed (i.e. com- 
pact without boundary). In the framework of the 3+1 
formalism of general relativity, one may think of each 
surface St as being the intersection of TL with a spacelike 
hypersurface S t arising from some 3+1 foliation (E t ) te K 
of M.: S t — TL n S t . Such a viewpoint will be called 
hereafter a 3+1 perspective (e.g. [12]). It will not be 
used in the mainstream of this article, except for making 
remarks and connections with previous works. Indeed, 
we will deal only with quantities intrinsic to TL and its 
foliation (S t ) te R. Besides, let us recall that for a dynam- 
ical horizon the foliation (St) tern, by marginally trapped 
surfaces is unique (up to a relabeling 1 1— > t') [19]. 

Demanding that the 2-surface St is spacelike amounts 
to saying that the metric q induced by the spacetime 
metric g onto St is positive definite (i.e. Riemannian). 
In particular q is not degenerate and at each point p € St, 



the following orthogonal decomposition holds: 

T P (M) =T p (S t )®T p (S t )- L , (2.1) 

where T P (M) [resp. T p (S t )] denotes the space of vec- 
tors tangent to M [resp. to St] at the point p, and 
T^St) 1 ^ denotes the space of vectors orthogonal to St 
at p. Both vector spaces T p (St) and T p (St) ± are two- 
dimensional. Let us then denote by q the orthogonal 
projector onto T p (S t ): q(v) = v v G T p (S t ) and 

q(v) = v € T p (St) ■ In this article, we shall take 

a 4-dimensional point of view on the induced metric q by 
setting q(u,v) — if any of the two vectors u and v in 
T p (M) belongs to T p (S t ) ± - Then, if in a given basis, the 
components of q are q a p, the components of q are q a p, 
where the index a has been raised with the metric g. 

Given a generic tensor A on T p (A4) of covariance type 
(m, n), we define a new tensor q* A of the same covariance 
type thanks to the projector q: 

A^-^ Ux .. Mn . (2.2) 

Note that for a vector q*v = q(v) and for a 1-form, 
q*u> = u> oq. Note also that for multilinear forms intrin- 
sic to the 2-dimensional manifold St, q* can be viewed 
as the "push-forward" operator which transforms them 
to multilinear forms acting on the 4-dimcnsional space 
T p (A4) (for vectors and more generally contravariant ten- 
sors, the push-forward operator is canonically provided 
by the embedding of St in M). A tensor A on M will 
be said tangent to St if q*A = A. 

B. Evolution vector 



1 Latin indices from the beginning of the alphabet (a, b, ...) run 
in {2, 3}, Latin indices starting from the letter i run in {1, 2, 3}, 
whereas Greek indices run in {0, 1, 2, 3}. 



Let us denote by h the vector field on TL such that (see 
Fig. 1) (i) h is tangent to TL, (ii) at any point in TL, h 
is orthogonal to the surface St going through this point, 
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(iii) the length of h is associated with the parameter t 
labeling the surfaces (<S t ) by 



C h t=l, 



(2.3) 



where Ch denotes the Lie derivative along h. In the 
present case (scalar field t), we have of course Cht = 
h^dpX = (dt,h), where brackets are used to denote the 
action of 1-forms on vectors. Given the foliation (<S t ) tg R, 
the conditions (i), (ii) and (iii) define h uniquely. Note 
however, that if the leaves St are relabeled by a new 
parameter t' = F(t) (where F : M — > K is a smooth 
one-to-one map), then h is transformed into 



ti = [F'(t)]- 1 ^ 



(2.4) 



so that Ch' f = 1 . 

An immediate consequence of Eq. (2.3) is that the 2- 
surfaces S t are Lie-dragged by the vector field h: given 
an infinitesimal parameter St, the image of the surface S t 
by the displacement of each of its points by the vector 
5th is the surface St+st (cf. Fig. 1). For this reason, h 
is the natural vector field to describe the "evolution" of 
quantities across the foliation of H. In particular, we will 
consider the Lie derivative along h, Ch as the "evolution 
operator" 2 along H. Since h Lie-drags the surfaces (<S t ), 
it transports any vector tangent to St to a vector tangent 
to St+st- In other words, 



VveT(St), C h veT(S t ), 



(2.5) 



where T{S t ) denotes the space of vector fields defined 
on H and which are tangent to St- Although the vector 
field h is not tangent to St, we can use property (2.5) to 
extend the definition of Ch to 1-forms u> acting in T{St) 
(i.e. 2-dimensional 1-forms associated with the manifold 
structure of <S t ), by setting 

\/veT(S t ), (£ h u,v):=£ h (u,v)-(w,£hv). (2.6) 

Note that the right-hand side of this equation is well 
defined thanks to Eq. (2.5). The definition of Ch is then 
extended immediately to any tensor field on St via tensor 
products and Leibniz' rule, e.g. Ch (u>i(g)u>2) := Ch Wi® 
u>2 + Wi (g> £(,«2- Given a multilinear form field A on 
St , we then denote by s Ch A the push- forward (via the 
projector q) of the derivative Ch A defined above: 



s C h A:=q*C h A. 



(2.7) 



One can then show that (see Appendix A of Ref. [12] for 
details): 



s C h A = q*C h q*A, 



(2.8) 



2 The term "evolution" stands for "variation as t increases" and 
can be made more concrete is one adopts the 3+1 perspective 
mentioned in Sec. II A. 



where the Lie derivative in the right-hand side is the stan- 
dard Lie derivative along h within the manifold M.. 

Owing to the fundamental property (2.3) and the re- 
sulting Lie-dragging of the surfaces (S t ), it is not sur- 
prising that the vector h has been introduced by many 
authors when studying foliation of hypcrsurfaces, in vari- 
ous contexts: h and Ch were denoted respectively £ and 
D/dt by Damour [4-6] in his black hole mechanics (H 
was then taken to be an event horizon); h was denoted 

d\ n ^ by Eardley [22] in his study of black hole boundary 
conditions for 3+1 numerical relativity, since it was then 
viewed as the part of the evolution vector d/dt which is 
normal to St in a coordinate system adapted to Ji. Sim- 
ilarly, h is denoted Q by Cook [23] when searching for 
boundary conditions for initial data representing quasi- 
stationary black holes. More recently, in the context of 
trapping and dynamical horizons, h has been denoted V 
by Ashtekar and Krishnan [16], V by Booth and Fairhurst 
[24, 25] and £ by Hayward [17, 18]. 

Let C be the scalar field defined on H as half the scalar 
square of h: 



C:=-hh, 
2 



(2.9) 



where a dot is used to denote the scalar product taken 
with the metric g. Since h is normal to St, an orthogonal 
vector basis of T p (Ti) is {h,e2,e^), where {e2,e^) is an 
orthonormal basis of T p (St). In this basis, the matrix of 
the metric induced by g on H is diag(2C, 1,1). We then 
conclude that 

Ti is spacelikc -^=> C > <=+ h is spacelikc 

H is null C = h is null 

H is timelike C < h is timelike. 

(2.10) 



III. EXTRINSIC GEOMETRY OF A 
SPACELIKE 2-SURFACE 

In this section, we review some basic results about the 
extrinsic geometry of a single spacelike 2-surface — not 
necessarily a member of a foliation. For future purpose, 
we take care to provide rather general definitions, for in- 
stance not limiting the definition of expansion and shear 
to null vectors, as usually done, nor limiting the defini- 
tion of the normal fundamental forms to some privileged 
normal frame. 



A. Expansion and shear along normal vectors 

Let us consider a fixed 2-surface S t . We denote by 
T(S t ) ± the space of vector fields v which are defined 
on St and everywhere normal to St'- Vp G 5(, v{p) € 
T p (St) . For any v G T(5i) ± , we define the deformation 
tensor of S t along v as the bilinear form 



or 6 



(„) 

a/3 



= V 



feci p 



(3.1) 
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where V is the afhne connection associated with the 
spacetime metric g and the underlining is used to denote 
in an index-free way the 1-form v canonically associated 
to the vector field v by the metric g. Note that thanks 
to the projector q in Eq. (3.1), &( v ' is independent of 
the values of v away from St (some extension of v in an 
open neighborhood of St being required for the spacetime 
covariant derivative V« to be well defined) . It is easy to 
see that the bilinear form ©W is symmetric, as the con- 
sequence of v being normal to the surface S t (Weingarten 
property). 

Let us consider the metric q induced by g on the 2- 
surfaces deduced from S t by Lie-dragging along v (recall 
that q is a priori defined only on St] we have of course 

q = q). Taking into account the symmetry of ©^ 
and expressing the Lie derivative in terms of V yields 
q*C v q = q*V v q + 2® ( - v \ Now, from the idempotent 
character of q, it is easy to see that q*W v q — 0, so that 
finally one ends with 



1 



q*C v q. 



(3.2) 



This equality justifies the name deformation tensor given 
to ©("); ©M measures the variation of the metric in 
S t when this surface is Lie-dragged along the vector v. 
Decomposing ©W into a trace part and a traceless part 
results in the definition of the expansion rate of St along 
v: 



0<«) :=qf"Q$ = C v In VI, 
and the shear tensor of St along v: 



,(*) 
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(3.3) 



(3.4) 



In Eq. (3.3) the second equality results from Eq. (3.2), 
q being the determinant of the components q a b with re- 
spect to a coordinate system x a = (a; 2 , a; 3 ) of the induced 
metric on the surface obtained from St by Lie-drag along 
v. Since Jq is related to the surface element s e of St by 
s e = y/q dx 2 A dx 3 , we see by considering coordinates x a 
constant along v field lines that 9^ is nothing but the 
relative rate of change of the area of a surface element 
Lie-dragged by v from S t : 



C v s e = 0W s e, 
hence the name expansion rate given to 6^"\ 



(3.5) 



B. Normal frames 

From the decomposition (2.1) and the spacelike char- 
acter of St, we see that the restriction of the metric g 
to the vector plane 71 (St) orthogonal to S t must be of 
signature (—,+). There are then two natural choices of 
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FIG. 2: Vector plane T p (St) normal to St at a given point 
p, with some orthonormal frame (n, s) and some null frame 
(£., k). The directions of £ and fe are uniquely defined as the 
intersections of T p (St) ± with the light cone emanating from 
p, whereas the directions of n and s can be changed by an 
arbitrary boost in a direction normal to St- 



pairs of vectors for generating this plane: (i) an orthonor- 
mal basis (n, s), i.e. a timelike vector n and a spacelike 
vector s satisfying 



n ■ n = — 1, s ■ s = 1, n ■ s = 0; 



(3.6) 



(ii) a pair of linearly independent future-directed null vec- 
tors {£., k); this choice is permissible since the signature 
(— , +) implies that T p (St) contains two null directions, 
which are actually the intersections of the null cone em- 
anating from p with T p (St) (cf. Fig. 2): 



£■£ = 0, kk = Q, Ik 



(3.7) 



where £ k is negative (hence written as minus some expo- 
nential) as a result of both £ and k being future-directed. 

In both cases, there is not a unique choice: in case (i), 
the timelike and spacelike directions can be changed by 
a boost in an arbitrary direction normal to St , leading to 
a new pair of basis vectors: 



n 
s' 



cosh rj n + sinh r\ s 
sinh rjn + cosh rj s, 



(3.8) 
(3.9) 



where r\ S M is the boost parameter. The choice (n, s) 
can be made unique by invoking some extra structure 
like the global foliation (St)tgg of M. arising from the 
3+1 perspective mentioned in Sec. II A, n being then 
the future directed unit normal to E t and s one of the 
two unit normals to St which are tangent to S{. An- 
other definite choice of (n,s) can be performed when Tt 
is spacelike (resp. timelike), by demanding that n (resp. 
s) is normal to 7i\ s (resp. n) is then colinear to the 
evolution vector h, and in particular lies in TL. This is 
the choice adopted by Ashtekar and Krishnan [10, 15, 16] 
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for dynamical horizons, which are always spacelike (see 
Appendix A). 

In case (ii), the two null directions are unique, but the 
vectors £ and k can be rescaled arbitrarily by 



e 
k' 



lik, 



A > 
fi > o, 



(3.10) 
(3.11) 



where the positive sign of A and /z is chosen to preserve 
the future orientation. One may reduce the arbitrariness 
by fixing the scalar product £ ■ k to —1 [choice a = 
in Eq. (3.7)], but this determines only /i as being A -1 
and leaves the degree of freedom on A. We will see in 
Sec. IV B that, when considering not a single St, but the 
whole foliation (<S t ) tgR , this ambiguity can be fixed in a 
natural way, leading to a unique choice of {£, k). 

Note that for the choice (i), the orthogonal projector 
q on S t is expressible as 



q = 1 + (n, .) n - (s, .) s 



(3.12) 



(equivalently q 
choice (ii) 



g + n® n — s®s), whereas for the 



q = l + e- a (k, .)£ + e~° ' (£, .) k 



(3.13) 



(equivalently q = g + e~ a k ® £ + e~ a £ ® k). 

St is called a trapped surface if, in addition of being 
spacelike and closed, it satisfies 6^ < and 9^ < 0, 
and a marginally trapped surface if 6^ = and 9^ < 0, 
or 9^1 < and 9^ = [26]. If one of the two null di- 
rections, £ say, can be selected as being "outgoing" , St is 
called an outer trapped surface if 9^ < (irrespectively 
of the sign of 9^) [13]. It is called a marginally outer 
trapped surface (MOTS) if 6»W = 0. Notice that all these 
definitions are unaffected by the rescaling (3.10)-(3.11) 
of the null vectors £ and k. 



Second fundamental tensor 



Let (n, s) be an orthonormal frame of T(5i) x ; inserting 
expression (3.12) for q a v in the above relation and making 
use of definition (3.1) leads to 



K~ a — r, a p)( n ' 

K, 01 - n u /?7 



/3 7 ■ 



(3.17) 



Similarly, if one uses instead a null frame (£, k) for 
T(S t ) ± , expression (3.13) for q a v leads to 



0i 



= e 



e 



)■ 



(3.18) 



It is clear on formulas (3.17) and (3.18) that the sec- 
ond fundamental tensor is orthogonal to St in its first 
index, and symmetric and tangent to St in its second 
and third indices. The reader more familiar with the hy- 
persurface case should note that the second fundamental 
tensor for a hypersurface writes, instead of Eq. (3.17), 



K, 1 



0i 



-n a Kf3 7 , where n a is the normal to the hyper- 



surface and Kp 7 its second fundamental form or extrinsic 
curvature tensor. 



D. Normal fundamental forms 

Contrary to the case of a hypersurface, the extrinsic 
geometry of the 2-surface S t is not entirely specified by 
the second fundamental tensor K. Indeed, because it 
involves only the deformation tensors ©(•) of the normals 
to S t [cf. Eqs. (3.17) and (3.18)], K encodes only the part 
of the variation of St's normals which is parallel to St- 
It does not encode the variation of the two normals with 
respect to each other. The latter is devoted to the normal 
fundamental forms which are the 1-forms defined by (cf. 
e.g. [14] or [29]) 



Cl {n) := s-Vff 



n 



or n<») := aM V v n"g"J(3.19) 



tt« := n ■ V q -s [or Q« := n„V v ^ <f a ] (3.20) 



As for any non-null submanifold of (Ai,g), the second 
fundamental tensor of S t (also called extrinsic imbedding 
curvature tensor [27] or shape tensor [28]) is defined as 
the tensor K, of type (1, 2) relating the covariant deriva- 
tive of a vector tangent to St taken with the spacetime 
connection V to that taken with the connection in St 
compatible with the induced metric q, hereafter denoted 
byX>: 

V{u,v) eT(S t f, V u v = T> u v + K(u,v). (3.14) 
From the fundamental relation 



T»A = q*VA, 



(3.15) 



valid for any tensorial field A tangent to St , it is easy to 
express 7C in terms of the derivative of q: 



(3.16) 



if one considers an orthonormal frame (n,s) of T(S t ) ± 
and by 



V~£ kV ^ 



or 



or 



i 



k p £p 



k»VJ»q" 



(3.21) 



1 



k p £p 



l^ v k»q\ 



(3.22) 

if a null frame (£,k) of T(St) ± is considered instead. 
Note that, thanks to the projector q, the definition of 
the normal fundamental forms does not depend upon the 
values of the normal fields n, s, £ and k away from the 
2-surface S t - Note also that thanks to the division by 
k ■ £, the value of Sl^ does not depend on the choice of 
the null vector k complementary to £. From the orthog- 
onality relations (3.6) and (3.7), we have the immediate 
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properties: 

n (a) = _ n (n) (3 23 ^ 

n (fc) = _ n W +5Jffi (3.24) 

We can also relate the (•£, fe)-type normal fundamental 
forms to the (n, s)-type ones by choosing the canonical 
null frame associated with a given orthonormal frame 
(n, s), namely £ = n + s and k = n — s. Then 

n W = (Tl) [^ = n + s] (3.25) 
n (fe) = _ n (») [fc = n _a]. (3.26) 

Note that, if one considers a non-null hypersurface in- 
stead of a 2-surface, the analog of definition (3.19) would 
be n( n ) := n ■ V^n, since there is only one normal n. 
But this expression vanishes identically by virtue of the 
normalization of n (n ■ n = 1 for a timelike hypersur- 
face, and —1 for a spacelike one). Consequently, the ex- 
trinsic curvature of a non-null hypersurface is entirely 
characterized by the second fundamental form K. For a 
null hypersurface, with normal £, the orthogonal projec- 
tor q is not defined (as a result of q being degenerate) . 
The relevant quantity is then the 1-form fiO defined 
by Eq. (3.21) but with q substituted with the orthog- 
onal projector to some spacelike 2-surface embedded in 
the hypersurface and k substituted with a transverse null 
vector. It is then called the Hdjicek 1-form [30, 31] (see 
also [12]). 

The normal fundamental forms can be interpreted in 
terms of the connection 1-forms associated with respect 
to some tetrad. Indeed let e a = (n, s,e 2 ,e 3 ) be an or- 
thonormal tetrad [(e2,e3) is then an orthonormal basis 
of T p (St)\. The connection 1-forms associated with this 
tetrad are the 1-forms u}@ a such that for any vector field 
v on A4, V B e a = (u>» a ,v) e p . Then, from Eq. (3.19), 

fiW = q*u\. (3.27) 

An equivalent phrasing of this is saying that the two non- 
trivial components of fi( n ) with respect to the dual frame 
(e Q ), namely (a = 2, 3), are identical to some of the 
connection coefficients T a p associated with the tetrad 
(e Q ): 

= TV (3.28) 

Relations (3.27) and (3.28) justify the alternative names 
external rotation coefficients [27] and connection on the 
normal bundle [25, 29, 36] given to the normal funda- 
mental forms. 

The normal fundamental forms depend on the normal 
frame. Indeed a change of normal frame (n, s) i— > («.', s') 
according to Eqs. (3.8)-(3.9) leads to 

n (n "> =n (n) +Vn, (3.29) 

whereas a change of null normal frame (£,k) i— ► {£',k') 
according to Eqs. (3.10)-(3.11) leads to 



On the contrary, the second fundamental tensor K, in- 
troduced in the previous section does not depend on 
the choice of the normal frame: this is obvious from 
Eq. (3.16) which involves only the projector q, and this 
can be checked easily from the expressions in term of the 
normal frames [Eqs. (3.17) and (3.18)], by substituting 
the transformation laws (3.8)-(3.9) and (3.10)-(3.11). We 
refer the reader to Carter's article [27] for an extended 
discussion of this dependence of the normal fundamental 
forms with respect to the normal frames. 

IV. EXTRINSIC GEOMETRY OF THE 
FOLIATION 

Section III introduced quantities relative to a single 2- 
surface S t . Here we investigate quantities defined with 
respect to the family (S t )tm foliating Ti. 

A. Dual-null description of the foliation (St)tem 

A very convenient way to study the foliation (St)tev. 
is to employ of the dual-null formalism of Hayward 
[14, 32, 33] (see also [34, 35]), which we recall here, 
adapting the notations to our purpose (see Table I for 
the correspondence with Hayward's notations). 

Let us consider, in the neighborhood of TL, two fam- 
ilies of null hypersurfaces, {U u ) ue m. and (V„)„ e R, which 
intersect in spatial 2-surfaces such that each 2-surface S t 
is one of these intersections, i.e. for each ! e 1, there 
exists a value of u, uo(t) say, and a value of v, v (t) say, 
such that St is the intersection between U uo i t ) and V vo i t ) : 

VteE, 5 t =W Bo(() nv„ o(t) . (4.1) 

Such a dual- null foliation always exists, hi Ua (t) (resp. 
V„ (t)) being nothing but the hypersurface generated by 
light rays outgoing (resp. ingoing) orthogonally from S t . 
Moreover, if TL is spacelike or timelikc, the dual-null foli- 
ation is unique. If TL is null, then TL coincides with U Uo m 
and u (t) = const.; there is then the degree of freedom 
of choosing the foliation (U u ) outside of TL. 

Let £ and k be the null normal vectors to respectively 
U u and V v and dual (up to a sign) to the gradient 1-forms 
du and dv: 

£:=-du and fe := -dv. (4.2) 

Since St belongs to both U Uo tt) and V„ ( t ), both £ and k 
are normal to St, and therefore constitute a null frame 
of T(S t ) ± , similar to those considered in Sec. IIIB. Let 
us denote f_ the scalar field o associated to the scalar 
product of £ and k by Eq. (3.7): 

£-k=:~e s . (4.3) 

From the definition (4.2), the 1-forms £ and k are closed: 



=nW+X>lnA. 



(3.30) 



d£ = and dfc = 0. (4.4) 
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The vectors £ and k being null, it follows immediately 
from Eq. (4.4) that 



V g £ = 



and 



0. 



(4.5) 



i.e. the field lines of £ and k are geodesies: they are 
the light rays emanating orthogonally from S t . Hayward 
[14, 17, 18, 32, 33] introduces another pair of null vectors 
by setting 



£ :-- 



-n 



and 



k :-- 



- f k. 



(4.6) 



These vectors have the fundamental property of Lie- 
dragging the hypersurfaces V v and U u respectively, i.e. 
they obey 



1 



and 



1. 



(4.7) 



but contrary to £ and k, the 1-forms £ and k are not 
closed, since we deduce from Eq. (4.4) that 



d£-- 
Besides, 



-d/A£ 



and 



dfc = -df Afe. 



£k = -e~ f . 



(4i 



(4.9) 



The anholonomicity 1-form, also called twist 1-form, 
is defined by [cf. Appendix B of Ref. [14], in conjunction 
with our definitions (3.21) and (3.22)]: 



■BE7 



: = 1 (n* - nw) . (4.io) 

From Eq. (3.24) with a = /, we get 

vj = -nw + ix>/ = n< s > - ix>/. (4.ii) 

According to the scaling law (3.30) with A = [cf. 
Eq. (4.6)], we can re-express the anholonomicity 1-form 
in terms of the normal fundamental forms associated with 
£ and k : 



- ix>/ = + ^Vf. (4.12) 



Thanks to Eq. (4.4), we can easily re-express vj in term 
of the commutator of k and £, as well as that of k and £: 



-f 



of 



■uj = 



q [k,£] = -q [k,£] 



—f f 
or w a = —q a ^,£Y = —q a ^[k,lY 



(4.13) 



(4.14) 



This justifies the terms anholonomicity and twist given to 
vj: according to Frobenius theorem, the 2-planes T p (St) 
are integrable in 2-surfaces when t varies, if, and only if, 
the commutator of two generating vectors, e.g. k and £, 
satisfies [k,£] € T p (S t ) ± ; from Eq. (4.13), this is equiva- 
lent to vj = 0. 



TABLE I: Correspondence between our notations and those 



of Hayward. 


Llllb WU1K. 


IT d. y W a,L CI l-L^J 




Q 


ft 


ft 






7-1 


u 




a; - 


V 


s+ 


X 


2 


iV — 




k 


A+ 


1 / _L \ 

9 1 (n+) 


£ 


n_ 




k 


n+ 




f 


/ 


/ 


£ 


— r+ = e ^Af_ 


7 

«+ 


k 


M- — r- — e 3 N+ 


Z- 




LU 


LU 




— p-f ft 


M+) 


n (fc) 


-e- / /3+ 


C(-) 


n<*> 




-<(-) 


n<*> 




-<(+) 


h 






B 






C/A 




-r 


m 




r 



B. Normal null frame associated with the evolution 
vector h 

The evolution vector h introduced in Sec. II B belongs 
to the plane orthogonal to S t - Following Booth and 
Fairhurst [24, 25], we notice that there exists a unique 
pair of null vectors (£, k) in that plane such that (sec 
Fig. 3) 



h = £-Ck 



and 



£k = -1, 



(4.15) 



where C is related to the scalar square of h by Eq. (2.9). 
Thus we may say that the foliation (St) tern, entirely fixes, 
via its evolution vector h, the ambiguities in the choice 
of the null normal frame (£, k) discussed in Sec. IIIB. 

The vectors £ and k are necessarily colinear to respec- 
tively the vectors £ and k associated with the dual-null 
foliation introduced above, i.e. there exists two positive 
scalar fields, A and B, such that 



£ = A£ and k = Bk. 



(4.16) 



Actually, we will use Eq. (4.16) to define £ and k away 
from TC, Eq. (4.15) defining them only on Ti. . However, all 
the results presented here are independent of the values 
of A and B away from Ti.. The normalization £ ■ k = — 1, 
combined with Eq. (4.3) relates the product AB to /: 

AB = e- f or / = -ln(AB). (4.17) 

Then, from Eqs. (4.16) and (4.6), 

£ = B~ x i and k = A' 1 ft, (4.18) 
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C. "Surface-gravity" 1-forms 



FIG. 3: Null vectors (£., fc) associated with the evolution vec- 
tor h by h — I — Ck; the plane of the figure is the plane 
T p (S t ) ± , so that S t is reduced to a point. 



which implies 



h = B-H- CA- X k. 



(4.19) 

Consequently, taking into account Eqs. (4.2), (4.6) and 
(4-7), 



C h u = -CA- 1 



and Ch v = B 1 . 



(4.20) 



On the other side, since u = uo(t) and v = vo(t), 

Ch u = w'o(0 and Ch v = v'o{t), where we have used the 
fundamental property Ch t = 1 defining h [Eq. (2.3)]. 
We then conclude that 



B U l/v' (t) 



and 



C/A=-u' (t). (4.21) 



This implies that, on H, the fields B and C/A are func- 
tions of t only; in particular, they are constant on each 
2-surface S t : 

T>B = and T>(C/A) = 0. (4.22) 
Using Eqs. (4.16) and (4.4), we get 
d£ = d\nAA£ and dfc = dlnB Afc, (4.23) 
from which we obtain 

V £ £ = v (t) £ with V( t y.= CebxA, (4.24) 



Vfc fc = V( k) k 



with 



v {k) :=C k hiB. (4.25) 



and Vt k \ are the inaffinity parameters of the null 
vector fields £ and k. Using the definitions (3.1), (3.21) 
and (3.22), we then get an expression for the spacetime 
gradients of £ and fc: 

VI = © w + I <g> fl w - v w l <g> fc - V k l <g> £ (4.26) 

Vfe = (fe) - fc ® fiW - i/ (fe) fe <g> £ - V t k <g> fc, (4.27) 

where we have used Eq. (3.24) with a = to set fl^ = 
-fiW. Besides, from Eqs. (4.23) and (4.22), we get useful 
identities: 



<fV fe £ = -fJ w +X>lnA 



(4.28) 
(4.29) 



Let us define the 1-form 
1 



fc 



kJ^^ q 



(4.30) 



where q denotes the orthogonal projector on the vector 

_L 

plane T p {St) , i.e. the complementary of q: 1 = q+ q. 
The definition (4.30) is similar to the definition (3.21) of 

_L 

S~2( £ ), except for q replaced by q. Hence, whereas II" 
was defined for a single 2-surface St, requires the 
knowledge of the null normal £ in directions normal to 
S t . From Eq. (4.24) and (4.25), the inaffinity parameters 
V((\ and are recovered by applying the 1-form tv® 
to respectively £ and — fc: 



{k^ ,£) = vii} and 
A useful relation is then 



,fc) 



(K [t, ,h) = v {t) + Cv {k) . 



(4.31) 
(4.32) 



D. Trapping horizons and dynamical horizons 

Let us recall here the various definitions involved in 
the local characterizations of black holes mentioned in 
the Introduction. The hypersurface Ti. equipped with the 
spacelike foliation (St)teR is called a marginally outer 
trapped tube (MOTT) [11] if each leaf St is a marginal 
outer trapped surface (cf. Sec. IIIB), i.e. if 6 = at any 
point in Ti. Following Hay ward [14] a trapping horizon is 
a MOTT on which 9^ ± and C k 9^ ^ 0, being quali- 
fied as future trapped horizon if 6^ k ' < and future outer 
trapped horizon (FOTH) if 9^ < and C k 9^ < 0, 
the latter subcase being the one relevant for black holes 
(see Ref. [11] for a discussion). The dynamical horizons 
introduced by Ashtekar and Krishnan [10, 15, 16] are 
MOTT such that (i) Ti is spacelike and (ii) 9^ < 0. In 
particular, a spacelike future trapping horizon is a dy- 
namical horizon. When it is null, a MOTT is called a 
non-expanding horizon [30, 31, 37]. It corresponds to a 
black hole in equilibrium. 



V. THE GENERALIZED 
DAMOUR-NAVIER-STOKES EQUATION 

A. Original Damour-Navier-Stokes equation 

In the case where the hypersurface Ti is null, and in 
particular when Ti is the event horizon of a black hole, 
h = £ and the Damour-Navier-Stokes equation [5, 6, 12, 
38] writes 

1 



s cM t] 



T>v ( , 



'(/) 

+8'Kq t T ■ £. 



(5.1) 
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FIG. 4: Same as Fig. 3 but with in addition the vector m 
normal to the hypersurface TL. 

This equation is derived from the Einstein equation, as 
the presence of the stress-energy tensor T testifies 3 . It 
has exactly the same structure as a 2-dimensional Navier- 
Stokes equation: dividing Eq. (5.1) by 87r, — H^/(8tt) 
is interpreted by Damour [5, 6] as a momentum surface 
density, iy^/(8n) as a "fluid" pressure, l/(167r) as a shear 
viscosity (er^ being the shear tensor), — l/(167r) in front 
of £>6>W as a bulk viscosity and q*T ■ t as a force surface 
density. The reader is referred to Chap. VI of Ref. [8] for 
an extended discussion of this "viscous fluid" viewpoint. 



B. Derivation of the generalized equation 

First of all, it must be noted that in the Damour- 
Navier-Stokes equation (5.1), the vector field t plays two 
different roles: it is both the evolution vector along TL 
(obviously in a term like s Ct ) and the normal to TL (in 
a term like q*T ■ £). When TL is no longer null, these two 
roles have to be taken by two different vectors. We have 
already seen that the privileged evolution vector along 

3 We are using geometrized units, in which both the speed of light 
c and the gravitation constant G are set to 1 



TL is the vector h associated with the foliation (S t ) te R. 
Regarding the normal vector, it is natural to consider 

m:=£ + Ck. (5.2) 

Indeed this vector is normal to TL, since by construction 
m E T p (St) ± and m ■ h = 0, and in the limit where TL 
is null, it reduces to t. It can be viewed as the unique 
normal vector to TL whose projection onto TL along the 
ingoing null direction k is h (see Fig. 4). Note that the 
scalar square of m is the negative of that of h: 

mm = -2C = -hh. (5.3) 

A generalization of the Damour-Navier-Stokes equa- 
tion to the non-null case should contain the term q*T m, 
instead of q*T-£, in the right-hand side. By virtue of Ein- 
stein equation and the fact that q(m) = 0, &Ttq*T ■ m — 
q*R- m, where R is the Ricci tensor associated with the 
spacetime metric g. The starting point for getting the 
generalized Damour-Navier-Stokes equation will be then 
the contracted Ricci identity applied to the vector m and 
projected onto St- 

(V^m" - V„V,m") q\ = R^m»q\. (5.4) 

Now, by combining definition (5.2) with expressions 
(4.26) and (4.27), 

-IJfV^ - Ck a £ a V a kP + V Q C*fc /3 .(5.5) 

Substituting Eq. (5.5) for Vj,m M and V^m^ in Eq. (5.4), 
expanding and making use of identities (4.28), (4.29) and 
(4.32) yields, after some rearrangements, 



R^ m » q \ = <f Q v M e("^ + cfjfVprtp + ej,"'" n<f> + n£) - v a [e^ + («w , h)j - e<f>p v„ in a 

-ej fc >" V^C + 9^V a C + v {k) (V a C - CV a In A) , (5.6) 

=o 

where the "= 0" results from Eq. (4.22). Now, from the relation (3.15) between the derivatives T> and V, one has 
[making use of identities (4.28) and (4.29)] 

q\V„G^ = V^Q(™> + S^V^lnA. (5.7) 

Besides, expressing the Lie derivative in terms of V gives 

q \c h nw = q \h^^ + e™» . (5.8) 

Thanks to Eqs. (5.7) and (5.8), Eq. (5.6) reduces to 

R pv m»q\ = x^e^)" + q \c h + - v a (V™> + («w h)) + e^v a c, (5.9) 
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where Eq. (4.22) has been used to set to zero the term 0^ (CV^ In A — D^C) which had appeared. Now, from 
Eq. (2.8) and the property g*fiW = , the term g*£ h «W which appears in the above equation is nothing but 
s Ch^l^\ Re-expressing in terms of the shear tensor er( m ) and the expansion scalar 0( m ) via Eq. (3.4) and 

taking account the Einstein equation then leads to 

s c h nw + = t>{kW,h) - d ■ a^ + -D9^ - e^vc + 8n^T ■ m. (5.10) 



This is the generalization of Damour-Navier-Stokes equa- 
tion to the case where the foliated hypersurface H is not 
necessarily null. In the null limit, C = 0,h = m = £ and 
we recover Damour's original version, i.e. Eq. (5.1). In 
the non-null case, it is worth to notice that the obtained 
equation is not much more complicated than Eq. (5.1): 
apart from substitutions of £ by either h or m, as dis- 
cussed above, it contains only one extra term: 9^DC. 

C. Change of normal fundamental form 

Let us rewrite the generalized Damour-Navier-Stokes 
equation in terms of the normal fundamental form f2^ ) 
associated with a generic null vector £', instead of £. Set- 
ting £' = X£, with A > 0, fl^ is related to «W by 



Eq. (3.30), from which we deduce 

s c h + n& = s c h + 9^ 

+D(C h lnX)+9 ( - h) DlnX, (5.11) 

where we have used s ChD In A = T>(Ch In A). Besides, 
we note that the 1-form k'*' transforms as follows 

= K W+VxInA, (5.12) 
q 

from which 

T>(K (e '\h) =V(K ( - e \h) +D(C h In A). (5.13) 

Combining Eqs. (5.10), (5.11) and (5.13), and using 
0(h) = e {t) _ CQ (k) yields 



s C h rt e "> + 0< fc ) Sl^ = D{K^'\h) - D ■ *( m ) + ^D9^ + 0^2? In A - 6»( fe ) (DC + CD In A) + 87r<f T • m. (5.14) 

If we chose A = A^ 1 , then by virtue of Eq. (4.22), the term DC + CXMnA vanishes identically Since A = A^ 1 
corresponds to £' = £ [cf. Eq. (4.16)], Eq. (5.14) reduces then to 

s C h + = D(n^\h) - D • <y( m > + ^D9 {m ^ - 9^D\nA + 8nq*T ■ m. (5.15) 

I 



In the case where H is not null (C ^ 0), another way 
to set to zero the term DC + CD\n A in Eq. (5.14) is to 
choose 



A = F(t)/|C|, (5.16) 



where F(t) is an arbitrary function of t, since then 
DC + CD\nX = CX>ln(A|C|) = CD\uF(t) = 0. The 
simplest choice F(t) = 1 corresponds to the following 
decomposition of the evolution vector: h = ±(C£' — k') 
(with ± corresponding to C > and C < respectively), 
to be contrasted with the decomposition (4.15). Actually 
this amounts simply to swapping the vectors £ and k. 



D. Application to trapping horizons 

If 7i is a trapping horizon (or more generally a MOTT, 
cf. Sec. IV D), then <?W = and Eq. (5.15) becomes 

s C h fl® + 0< h > fl® = D{ K &,h) - D • 

+ l-D9^ + 8tto*T ■ m. (5.17) 

This equation is structurally identical to the original 
Damour-Navier-Stokes equation [Eq. (5.1)]: apart from 
substitutions of £ by either h or m, it does not contain 
any extra term. The differences are that the original 
Damour-Navier-Stokes applies to a null Ti but with 9^ 
not necessarily zero, whereas Eq. (5.17) is valid for both 
H null or spacelike, but assumes - 0. 
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VI. ANGULAR MOMENTUM 

Traditionally the concept of angular momentum is a 
global one and requires the evaluation of a Komar inte- 
gral at spatial infinity, assuming M. to be asymptotically 
flat and endowed with an axisymmetric Killing vector (cf. 
e.g. Ref. [39]). However, by means of some Hamiltonian 
analysis, the concept of angular momentum can be made 
quasilocal, as a quantity associated with the interior of 
the spacetime region delimited by the hypersurface H. 
The prototype of such quasilocal formulation is Brown- 
York analysis [40] which will be taken as the starting 
point for our discussion. 

A. Brown- York angular momentum 

Let us assume that the hypersurface 7i is timelike and 
is axisymmetric, with the associated Killing vector tp ly- 
ing in the 2-surfaces S t - The definition of angular mo- 
mentum by Brown and York [40] is then 

J:= <f (3,?)% (6.1) 

where s e is the surface element of St associated with the 
induced metric q ( s e = ^fq dx 2 A dx 3 for any coordinate 
system x a = (x 2 ,x 3 ) on St, with q := dct q a b) and the 
momentum surface density 1-form j is defined as follows. 
Adopting a 3+1 perspective (cf. Sec. II A), let E t be 
a spacelike hypersurface intersecting TC in St- Denoting 
by 7 and K the induced metric and extrinsic curvature 
tensor of E t , j is expressible as 

ja = —^P^S^a, (6.2) 

where s is the unit spacelike normal to St which lies in 
E( and P is the momentum canonically conjugate to 7: 

p«/3 = y^:V7 {Kj af3 - K af3 ) . (6.3) 

Since K is related to the gradient of the timelike unit nor- 
mal to E t , n, by K a0 = -W^n v 7^7^, we get, by insert- 
ing Eq. (6.3) in Eq. (6.2) and comparing with Eq. (3.19), 

j = -^ n) - (6.4) 

By considering the null vector £' := n + s and combining 
the transformation laws (3.25) and (3.30), one gets 

j = -^{n^+V\n\), (6.5) 

where A is the scale factor relating £ to n + s: n + s = 
\£. Now, the Killing equation for the vector tp and the 
fact that tp e T(S t ) imply "D ■ tp = 0. Consequently 
tp-T> In A = 2?- (In A (p) is a perfect divergence, the integral 



of which over the closed surface St vanishes. Therefore 
substituting Eq. (6.5) for j into Eq. (6.1) yields 

j = -i|(nm, y )%. (6.6) 

This expression is in perfect agreement with the inter- 
pretation of —Sl^/(8n) as a momentum surface density 
performed in Sec. V. 

B. Generalized angular momentum 

It may be noticed that the timelike character of the hy- 
persurface Tt, assumed in Brown- York Hamiltonian anal- 
ysis [40], does not play any role in the expression (6.6) 
of the angular momentum. Actually the definition of an- 
gular momentum, based on Eq. (6.6), has been extended 
to null hypersurfaces by Booth [41] (in a generalization 
of Brown- York analysis) and Ashtekar et al. [42] (in the 
framework of isolated horizons). 

It is also worth to notice that the independence of 
the integral defining J with respect to the normal fun- 
damental form (i.e. f|( n ) or fl^) stems only from 
the divergence-free property of the vector tp, which is 
a condition weaker than that of being a Killing vector. 
Therefore, one may relax the later and follow Booth and 
Fair hurst's recent analysis [25] to introduce a generalized 
angular momentum as follows. Let us assume that the 
2-surfaces S t have the topology of § 2 . Let tp be a vector 
field in T(St) which (i) has closed orbits and (ii) has van- 
ishing divergence with respect to the induced connection 
V: 

T>tp = 0. (6.7) 

The angular momentum associated with tp is then defined 
by [25] 

%):=-i|(nW^) s e, (6.8) 

which is a formula structurally identical to formula (6.6). 
The main difference is that tp is no longer the Killing 
vector reflecting the axisymmetry of St and uniquely de- 
fined by the normalization of the orbit lengths to 2ir, but 
merely a divergence-free vector field. Consequently, J de- 
pends on the choice of tp. However formula (6.8) shares 
with formula (6.6) the independence with respect to the 
choice of the normal fundamental form thanks 
to the divergence-free character of tp. Indeed, under 
a change of null normal £' = \£, fl^ is changed to 
= n W + x>i n A [Eq. (3.30)] and since T> ■ tp = 0, 
tp ■ T> In A = T> ■ (In A tp) is a perfect divergence, the inte- 
gral of which on St vanishes. 

As a further justification of definition (6.8), it is shown 
in Appendix A that, when TL is a dynamical horizon, 
J(tp) agrees with the generalized angular momentum de- 
fined by Ashtekar and Krishnan [10, 16]. 
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C. Angular momentum flux law 

For any 1-form u) and vector field tp defined on H and 
both tangent to S t for all t £ M., the following identity 
holds: 



s £ h {(u,,tp) s e] 
( s £ h u + 0Wu,<p) s e 



(6.9) 



where the first equality results from the Lie transport of 
the 2-surfaces St by the vector field h (cf. Sec. II B) and 
in the second equality the relation Ch s £ = 0^ s e has 
been used [cf. Eq. (3.5)]. 



Applying the above identity to the 1-form u> = Sl^ 
and employing the generalized Damour-Navier-Stokes 
equation (5.10) leads to an evolution equation for the gen- 
eralized angular momentum J(tp) defined by Eq. (6.8): 



dt 



1 

ToTT 



s t 



a^: £ v q-29^tp-T>C 



e — 



8tt 



(6.10) 



The notation ':' stands for a complete contraction, whereas the double arrow means that the two indices of o^ 1 ™) have 
been raised with the metric g: : C v q = a ( - m ' lab X v q a b- The integrals involving the pure gradients T>{n^,h) 

and T>9 i - mS> have been set to zero thanks to the property T> ■ tp = 0. Besides, we have written tp ■ (T> ■ er^™ 1 )) = 
2? . (a( m ) ■ ipj — <x( m ) : Dip, with the integral of the divergence D ■ (<x( m ) • tp) being zero since St is a closed manifold 

and, thanks to the symmetry of the shear tensor er( m ), 2v (m ^> : T>p = a^ ah {V a <p h + V b ip a ) =<r( m ) : C v q. 

The last integral in Eq. (6.10) occurs to take into account a possible variation of p along the evolution vector h. 
To make the variation of J more meaningful, it is natural to demand that the vector field tp is transported by h: 



C h tp = 0. 



(6.11) 



From now on, we assume that tp obeys to both conditions (6.7) and (6.11). Note that if tp is a symmetry generator 
of H which is tangent to S t , these two conditions are satisfied (in particular Ch tp = —C v h = 0). Then Eq. (6.10) 
simplifies to 



sT(m , 9) . e __L£ 



: C v q - 20^ tp ■ V C 



e. 



(6.12) 



r 



In the case where H is a null hypersurface, then C 
and Eq. (6.12) reduces to 



T(£,tp)' 



q"e. 



1 

is, u >~ J s t 

(6.13) 

We recover here Eq. (6.134) of the Membrane Paradigm 
book [8], where the first term is interpreted as the vari- 
ation of J due to the flux of angular momentum carried 
by matter and electromagnetic field at H and the second 
term accounts for the shear viscosity of H. 

Let us consider now the case where 7i is either timelike 
or spacelike: C ^ and we may express 6K fc ) in Eq. (6.12) 
as 6>( fc ) = (6>W - 0W) /C [cf. Eq. (4.15)], so that 

-d^tp-VC = -eWp-T>\n\C\+9 {h) p-T>\n\C\. (6.14) 



Now, the properties (6.7) and (6.11) fulfilled by p imply 
that the vector field O^ip is divergence- free on (S t ,q): 



t>- (e^pj = o. 

This follows from the identity 



(6.15) 



C h V ■ tp + O^V ■tp = 'D-(c h tp + 0W(p) , (6.16) 



which can be easily established, for instance by con- 
sidering a coordinate system (t 7 x 2 ,x 3 ) on Ti such that 
h = d/dt. As a consequence of Eq. (6.15), the integral 
over St of O^tp ■ X>ln |C| vanishes. Taking into account 
Eq. (6.14), we deduce then that Eq. (6.12) can be written 



dt 



J(tp) 



s t 



T(m, tp) s e 



16tt 



s t 



a (m ^. C^q-26^tp-'D\n\C\ 



(6.17) 
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We deduce immediately from this expression that if TL is 
a timelike or spacelike MOTT, i.e. if 9^ = 0, the angular 
momentum variation law takes a very simple form: 

(6-18) 

In particular, the above formula holds for a spacelikc fu- 
ture outer trapping horizon 4 and for a dynamical horizon. 
We have established this relation by assuming C / 0. 
Now the relation obtained in the null case, Eq. (6.13), 
is identical to Eq. (6.18) since m = £ when C = 0. 
Therefore, collecting the two results, we conclude that 
Eq. (6.18) holds for a future outer trapping horizon of 
any kind: dynamical horizon (C > 0) or non-expanding 
horizon (C = 0). 

It is worth to note that Eq. (6.18) has the same form as 
Eq. (6.13), which has been established for generic null hy- 
persurfaces, not necessarily non-expanding horizons. In 
the case where TL is a dynamical horizon, m is timclike 
and since <p is orthogonal to m, — T(m, <p) represents a 
momentum density along the spatial direction <p. There- 
fore we can attribute the first term in right-hand side of 
Eq. (6.18) to the flux of matter and electromagnetic an- 
gular momentum across St- Regarding the second term, 
it clearly vanishes if St is axisymmetric, with (p as a sym- 
metry generator (C v q = 0) . 



D. Relation to previous angular momentum laws 

Booth and Fairhurst [24] have obtained the angular 
momentum flux law (6.18) is the case where TL is a slowly 
evolving horizon. A slowly evolving horizon is a future 
outer trapping horizon which is close (in a sense made 
precise in Ref. [24]) to an isolated horizon. Eq. (6.18) is 
then obtained by an expansion to second order in e, where 
e ~ \J1C is the small parameter which measures the devi- 
ation from an isolated horizon 5 . At this order, note that 
fj(m) j g re pi a ced by <r^ in Booth and Fairhurst's version 
[their Eq. (14)]. Another difference with these Authors is 
that they did not assume that <p is divergence-free, but 
that it is close to a Killing vector of (S t ,q). 

In the case where TL is a dynamical horizon, Ashtekar 
and Krishnan [16] have also derived an angular momen- 
tum balance law, but in a time-integrated form, so that 
it involves 3-dimcnsional integrals. Actually the relation 
derived in Ref. [16] does not assume that TL is a MOTT 
and is valid for any spacelike hypersurface. It turns out 
that if we integrate with respect to t the angular mo- 
mentum law (6.17), which is also valid for any spacelike 



4 If the null energy condition holds, a non-null future outer trap- 
ping horizon is necessarily spacelike [14]. 

5 Note that the definition of angular momentum in Ref. [24] has 
the opposite sign than ours. 



TL, we recover exactly Ashtekar and Krishnan's version. 
This is shown in Appendix A. 



VII. CONCLUSION 

We have established, by means of Einstein equation, an 
identity valid for any hypersurface TL foliated by space- 
like 2-surfaces (S t )tew- This equation has the same form, 
up to some additional term, as the 2-dimensional Navier- 
Stokes equation obtained by Damour [5, 6] for describing 
the dynamics of event horizons. The evolution vector, 
giving the time derivative of the effective momentum sur- 
face density is the vector h tangent to TL, orthogonal to 
the leaves S t and which transports them into each other 
(Lie-dragging) . The role of the momentum surface den- 
sity is played by the normal fundamental form il") of St 
associated with the outgoing null normal £ whose projec- 
tion along the ingoing null direction is h. The pressure 
term involves both vectors £ and h, as it is the "surface 
gravity" 1-form of £ acting on h. The vector defining the 
shear and the expansion involved in the viscous terms is 
the vector m normal to TL and whose projection onto TL 
along the ingoing normal null direction is h. The vector 
m gives also the external force exercised by matter and 
electromagnetic fields, if any. 

It must be noted that another outgoing null vector can 
be selected instead of £, such as a tangent £ to the hy- 
persurfaces of outgoing light rays emanating orthogonally 
from St, the dual 1-form of which is closed: d£ = 0. The 
key point is that all normal fundamental forms and "sur- 
face gravity" 1-forms differ only by a gradient and their 
interchange alters only slightly the generalized Damour- 
Navier-Stokes equation. 

If the hypersurface TL is null, the three vectors £, h 
and m coincide and the equation obtained here reduces 
to the original Damour-Navier-Stokes equation [5, 6]. If 
TL is a marginally outer trapped tube, and in particular 
if it is a dynamical horizon or a future outer trapping 
horizon, the obtained equation, written in terms of £, is 
as simple as the original Damour-Navier-Stokes equation, 
the generic additional term vanishing in this case. 

The generalized Damour-Navier-Stokes equation has 
been used to derive a balance law for the angular momen- 
tum associated with each of the leaves St and a generic 
divergence-free vector. When TL is a dynamical horizon, 
this law is a time differential form of the law obtained 
by Ashtekar and Krishnan [10, 16]. When TL is a slowly 
evolving horizon, we recover the angular momentum flux 
law obtained by Booth and Fairhurst [24]. 



Acknowledgments 

It is a pleasure to thank Jose Luis Jaramillo for numer- 
ous discussions and to acknowledge the warm hospitality 
of the Yukawa Institute for Theoretical Physics and the 



14 



TABLE II: Correspondence between our notations and those 
of Ashtekar and Krishnan 



this work 



7 
K 

<1 
T> 

0« 

£ 
k 

1 dR. 
s/2C dt 

1 dR. p _ A dR. 
C dt C C dt ' 



Ashtekar and Krishnan [16] 



T 

f 

<7 

K 

£> 

K 

k 



w 

C 



Department of Physics of Kyoto University where this 
work was completed. 



null normal frame (see Table II): 

£:=n + s and k:=n — s. 



(A2) 



From Eqs. (Al), (4.15) and (5.2), we get immediately the 
relation between these vectors and the null vectors (£, k) 
associated with h and introduced in Sec. IV B: 

and k = V2Ck. (A3) 
Ashtekar and Krishnan have introduced the 1-form 

= V M n„^ a = fiW, (A4) 

hence W is nothing but the normal fundamental form 
fi( ft ) [cf. Eq. (3.19)]. 

Another 1-form introduced by them is 6 



(A5) 



Replacing s and £ by their respective expressions (Al) 
and (A3), and making use of Eq. (4.15), yields 



APPENDIX A: LINK WITH THE 3+1 
DESCRIPTION OF DYNAMICAL HORIZONS 

1. Basic relations 

When H is a dynamical horizon, it is a spacelike hyper- 
surface and Ashtekar and Krishnan [10, 16] have made 
use of the standard 3+1 formalism to describe it, by in- 
troducing its unit timelike future directed normal n, its 
positive definite induced metric 7 and its extrinsic curva- 
ture tensor K. Here we have put bars on the symbols de- 
noting them to stress that these objects are relative to H 
itself and not to some spacelike hypersurface S intersect- 
ing H in a 2-surface (as in the 3+1 perspective mentioned 
in Sec. II A). Note that our sign convention for the ex- 
trinsic curvature is the opposite of that of Ashtekar and 
Krishnan and that we are using n for the unit normal 
denoted f by Ashtekar and Krishnan. We also denote 
by s the unit spacelike normal to St lying in TL. (n, s) 
constitutes then an orthonormal frame normal to St- s 
is denoted by r by Ashtekar and Krishnan, but we priv- 
ilege here notations consistent with those introduced in 
Sec. IIIB. The correspondence between both sets of no- 
tations is given in Table. II. 

The normal n is necessarily colinear to m. Similarly 
s is necessarily colinear to h. From the norm of m 
[Eq. (5.3)] and h [Eq. (2.9)], we deduce 



1 



n 



2C 



and 



1 



2C 



h. 



(Al) 



Let us recall that C > for a dynamical horizon. 
Ashtekar and Krishnan [10, 16] consider the following 



Ca = ^"V„4 q\ - fc"V„£, q\. (A6) 

Thanks to the identities (4.24) and (4.28), and to the fact 
that for C > 0, Eq. (4.22) implies T>\nA = X>lnC, we 

get 



C = n (£) -DlnC. 



(A7) 



Now, since £ = A£, we have from the scaling law (3.30), 
fiW = nW+VlnA = fiW+XUnC. Hence we conclude 
that the quantity £ introduced by Ashtekar and Krish- 
nan [10, 16] is nothing but the normal fundamental form 
associated with the null vector £: 



C = n®. 



(A8) 



In Ashtekar and Krishnan analysis [10, 16], a privi- 
leged role is played by the area radius R, i.e. the scalar 
field on 7i, which is constant on each 2-surface St and 
related to the area a of this surface by a — 4nR 2 . An 
associated quantity is the lapse Nr defined as the norm 
of the gradient of R within TL: 



N R 



fidiRdjR. 



(A9) 



R is a function of t and we may write diR = (dR/ dt) dit 



(dR/dt) SijylC. From the normalization Y Js i s j = 1 



6 actually they introduced it as a vector, but we consider it here 
as a 1-form via the standard metric duality 
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and the positivity of dR/dt (area increase law [16]), we 
then obtain 



N R 



1 dR 
~W~dt' 



(A10) 



The null evolution vector considered by Ashtckar and 
Krishnan is = N R £. From Eqs. (A3), (A10), and 
(4.16), we can re-express it as 

1 dR „ A dR - / & \ 

*™ = clt l= clt L (A11) 

Notice that thanks to the property (4.22), the coefficient 
A/C dR/dt in front of £ is constant on each 2-surface S t - 



From Eq. (A4) and q(<p) = <p, we deduce immediately 
that 



1 

"8^ 



(A13) 



If we suppose now that ip is divergence-free with respect 
to the connection in S t - *D ■ ip — 0, then, by means of 
the transformation laws of normal fundamental forms, 
Eqs. (3.25) and (3.30), it is easy to see that J(<p) co- 
incides with the generalized angular momentum J(ip) as 
defined by Eq. (6.8): 



2. Angular momentum 



T> ■ if = 



J(<p) = J(<p). 



(A14) 



Ashtekar and Krishnan [10, 16] define the generalized 
angular momentum associated with a section S t and a 
vector field ip tangent to St by 



J(<p) :- 



K{<p,a) 



(A12) 



Regarding the angular momentum flux law, Ashtekar 
and Krishnan [16] have derived an integrated version of it 
from the momentum constraint equation relative to the 
hypersurface H. It writes 



J(ip,t 2 ) - J(<p,h) 



ah lo7r 



f [IK?* - K**)C V %] H e, 

J AH 



(A15) 



where AH. is a portion of 7i delimited by two surfaces, St ± and St 2 say, and n e is the volume 3-form on Ti associated 
with the metric 7. Note that we have restored the explicit dependence of J(ip) on St by writing J((p, t). Note also that 
Eq. (A15) holds for any spacelike hypersurface H., not necessarily a dynamical horizon. Let us express the integrand 
in the second integral in the right-hand side of Eq. (A15) in terms of fields defined on the 2-surfaces St- First of all, 
performing an orthogonal 2+1 decomposition of K with respect to St yields 



(A16) 
(A17) 



Besides, Cpjij = Dupj + Djpi, where Di is the covariant derivative associated with the 3-metric 7 on H and 
Pi '■= i w i tn 

DiVj = Vifj - ®fk<f k s 3 + SiS k D kVj . (A18) 
Using the property Ch <p = [Eq. (6.11)] and the relation (Al) between h and s, we can rewrite the above expression 



as 



D m = Vnpj + - sj-ej?) <p k + \p k V k InCsiSj. 

From Eqs. (A16), (A17) and (A19) and the divergence-free property of <p [Eq. (6.7)], we get 

{Kf j - K ij )C v % = 2<j {n)ab V aVb - e {n) p a V a \nC = ^ (fi) : C v q - 9 (n) ip T>\nC. 



(A19) 



(A20) 



Now from Eq. (Al), cr^ = tr^/V^C and 6»( n ) = 6»( m )/\/2C, so that we can write Ashtekar and Krishnan's 
integrated flux law (A15) as 

■/(*>, t 2 )- J{<p,h) = - ( -l=T(m^)"£-i-/ -±=\*W:£ v q-6W<p-'DlnC] "e. (A21) 
J AH V 2C J-D7T j AH V2C L J 
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On the other side, if we integrate in time our flux law (6.17), which, as Eq. (A21), is valid for any spacelike hypersurface 
H, we get 

J(<p,t 2 ) - J(<p,h) = - [ T(m, tp) n dt A s e - -i- / \v (m) : C v q - 29^<p ■ T>\nd H dt A s e, (A22) 
J An 167r J ah 1 J 



where H dt denotes the gradient 1-form of the scalar field 
t within the manifold 7i. Besides, from the basic prop- 
erties Cht = 1, h orthogonal to St and h ■ h = 2C, 
we deduce easily that h := 7 • h = 2C n dt, hence 
I := 7 • s = h/V2C = \f2C H dt. Now, since s is the 
unit vector orthogonal to S t , n e = sA s e, so that we have 

n e = V2C H dt A s e. (A23) 

This relation shows the equivalence of Eqs. (A21) and 
(A22), except at first glance for the 9^ term in 



Eq. (A21) which is replaced by 26>W in Eq. (A22). How- 
ever, 9 (m *> = 26»W - 0< h ) [cf. Eqs. (4.15) and (5.2)] and 
the integral over S t of 6»( h V -T>\tiC vanishes since the 
vector field 9^<p is divergence-free [Eq. (6.15)]. This 
proves that Eq. (A22) is identical to Eq. (A21), i.e. that 
for a spacclike hypersurface, and in particular for a dy- 
namical horizon, the integrated version of our angular 
momentum flux law (6.17) results in Ashtekar and Kr- 
ishnan [16] angular momentum balance equation. 
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